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POSITIVE HARRIS RECURRENCE AND EXPONENTIAL 
ERGODICITY OF THE BASIC AFFINE JUMP-DIFFUSION 


PENG JIN, BARBARA RUDIGER, AND GHIRAZ TRABELSI 

Abstract. In this paper we find the transition densities of the basic affine 
jump-diffusion (BAJD), which is introduced by Duffie and Garleanu [D. Duffie 
and N. Garleanu, Risk and valuation of collateralized debt obligations, Finan¬ 
cial Analysts Journal 57(1) (2001), pp. 41-59] as an extension of the CIR 
model with jumps. We prove the positive Harris recurrence and exponen¬ 
tial ergodicity of the BAJD. Furthermore we prove that the unique invariant 
probability measure tt of the BAJD is absolutely continuous with respect to 
the Lebesgue measure and we also derive a closed form formula for the density 
function of tt. 


1. Introduction 


In this paper we study the basic affine jump-diffusion (shorted as BAJD), which 
is given as the unique strong solution X := {Xt)t>o to the following stochastic 
differential equation 


dXt=a{e-Xt)dt + a^/YtdWt + dJt, Xq > 0, 


( 1 . 1 ) 


where a, 6, a are positive constants, (Wt)t>o is a 1-dimensional Brownian motion 
and (Jt)t>o is an independent 1-dimensional pure-jump Levy process with the 
Levy measure 



cde '^'^dy, y >0, 

0, y <0, 


for some constants c > 0 and d > 0. We assume that all the above processes are 
defined on some filtered probability space (X)t>Q,P). 

The process X = (Xt)t>o given by (1.1) has been introduced by Duffie and 
Garleanu [3] to describe the dynamics of default intensity. It was also used in [4] 
and [14] as a short-rate model. Due to its simple structure, it is later referred as 
the basic affine jump-diffusion. The existence and uniqueness of strong solutions 
to the SDE (1.1) follow from the main results of [6]. We should remark that the 
BAJD process X = (Xt)t>o in (1-1) stays non-negative, given that Xq > 0. This 
fact can be shown rigorously with the help of comparison theorems for SDEs, for 
more details we refer the reader to [6]. 
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As its name implies, the BAJD belongs to the class of affine processes. Roughly 
speaking, affine processes are Markov processes for which the logarithm of the char¬ 
acteristic function of the process is affine with respect to the initial state. Affine 
processes on the canonical state space x R” have been thoroughly investigated 
by Duffie et al [2], as well as in [13]. In particular, it was shown in [2] (see also [13]) 
that any stochastic continuous affine process on x R” is a Feller process and 
a complete characterization of its generator has been derived. Results on affine 
processes with the state space R+ can also be found in [4]. 

Affine processes have found vast applications in mathematical finance, because 
of their complexity and computational tractability. As mentioned in [2], these 
applications include the affine term structure models of interest rates, affine sto¬ 
chastic volatility models, and many others. 

Recently the long-term behavior of affine processes with the state space R+ 
has been studied in [14] (see also [12]), motivated by some financial applications 
in affine term structure models of interest rates. In particular they have found 
some sufficient conditions such that the affine process converges weakly to a limit 
distribution. This limit distribution was later shown in [11] as the unique invariant 
probability measure of the process. Under further sharper assumptions it was 
even shown in [15] that the convergence of the law of the process to its invariance 
probability measure under the total variation norm is exponentially fast, which 
is called the exponential ergodicity in the literature. The method used in [15] to 
show the exponential ergodicity is based on some coupling techniques. 

In this paper we investigate the long-time behavior of the BAJD. More pre¬ 
cisely, as the first main result of this paper, we show that the BAJD is positive 
Harris recurrent. As a well-known fact, Harris recurrence implies the existence 
of (up to the multiplication by a positive constant) unique invariant measure. 
Therefore our result on the positive Harris recurrence of the BAJD provides an¬ 
other way of proving the existence and uniqueness of invariant measures for the 
BAJD. Another consequence of the positive Harris recurrence is the limit theorem 
for additive functionals (see e.g. [10, Theorem 20.21]), namely (1/t) f{Xs)ds 
converges almost surely to f{x)'K{dx) as t ^ oo for any / G Sf,(R+), where 
TT is the unique invariant probability measure for the BAJD. Some applications 
of Harris recurrence in statistics and calibrations of some financial models can be 
found in [1] and [9]. 

As a key step in proving the positive Harris recurrence, we will derive a closed 
formula for the transition densities of the BAJD, which seems also to be a new 
result. In particular this formula indicates that the law of the BAJD process 
at any time t > 0 is a convolution of a mixture of Gamma-distributions with a 
noncentral chi-square distribution. We should point out that this fact has already 
been discovered in [4] for the BAJD in some special cases, namely for the BAJD 
whose parameters (a, 9, a, c, d) satisfy c/(a — a^d/2) G Z . In general, it is not easy 
to identify the transition density functions of affine processes. However, as shown 
in [5], a general density approximation procedure can be carried out for certain 
affine processes and in particular for the BAJD. Different from [5] we seek in this 
paper a closed formula of the transition densities of the BAJD. 
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Finally we show the exponential ergodicity of the BAJD. We should indicate 
that the BAJD does not satisfy the assumptions required in [15] in order to get 
the exponential ergodicity. Our method is also a different one and is based on the 
existence of a Foster-Lyapunov. 

The remainder of this paper is organized as follows. In Section 2 we collect some 
key facts on the BAJD. In Section 3 we introduce the so-called Bessel distributions 
and some mixtures of Bessel-distributions. In Section 4 we derive a closed formula 
for the transition densities of the BAJD. In Section 5 we first study some continuity 
properties of the transition densities of the BAJD and then show its positive Harris 
recurrence. In Section 6 we show the exponential ergodicity of the BAJD. 


2. Preliminaries 


In this section we recall some key facts on the BAJD. As defined in (1.1), it is 
the unique strong solution X = {Xt)t>o to the SDE 


dXt = a{9 — Xt)dt + a\/J^dWt + dJt, Xq > 0. 


Throughout this paper we denote Ex{-) and Px{') as the expectation and proba¬ 
bility respectively given the initial condition Xq = x, with a; > 0 being a constant. 

By the affine structure of the BAJD process X, the characteristic function of 
Xt (given that Xq = x) is of the form 


Ex [e“^*] = exp u) -|- x’tp{t, u)) , it G := {u S C : Jin < 0}, (2.1) 


where the functions (/)(t,u) and solve the generalized Riccati equations 


dt(l){t,u) = E{tl;{t,u)), (j){0,u) = 0, 

= R{'il){t,u)), ’tp{0,u) = u, 


( 2 . 2 ) 


with 


F{u) = aOu + 


d — u 


u 


eC\{d} 



— au, u € <C. 


By solving the system (2.2) we get 



(2.3) 


and 



where A = a — a‘^d/2. 
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According to (2.1), (2.3) and (2.4), the characteristic function of Xt is given by 


E*le"**l = I 


(1 - i «(1 - e -'))-". -■)-- ) ^ 


2a * 

■exp I 


( .-gA"-.., ). it A 4 0, 


vi). 


(d-u) ■ 

if A = 0. 

(2.5) 

Obviously [exp(uAt)] is continuous in t > 0 and thus the BAJD process X is 
stochastically continuous. 

We should point out that if we allow the parameter c to be 0, then the stochastic 
differential equation (1.1) turns into 

dZt = a{6 — Zt)dt + a\/^dWt, Zq = x > 0. (2.6) 


To avoid confusions we have used Zt instead of Xt here. The unique solution Z := 
{Zt)t>o to (2.6) is the well-known Cox-lngersoll-Ross (shorted as CIR) process and 
it holds 


= (1 - ^w(l - e “*)) ^ • exp (^ — 


xue 


- - e““*) 


(2.7) 


In Section 4 we will find a distribution vt on R+ such that 



if A ^ 0, 
if A = 0. 


( 2 . 8 ) 


Then it follows from (2.5), (2.7) and (2.8) that the distribution of the BAJD is 
the convolution of the distribution of the CIR process and Vf In light of this 
observation we can thus identify the transition probabilities p{t, x, y) of the BAJD 
with 

p(.t,x,y)=[ f{t,x,y-z)i^t{dz), x,y>0,t>0, (2.9) 

where f{t,x, y) denotes the transition densities of the CIR process. 


Remark 2.1. For a different way of representing the distribution of Xt as a con¬ 
volution we refer the reader to [4]. In fact it was indicated in [4, Remark 4.8] that 
the distribution of any affine process on R._|_ can be represented as the convolution 
of two distributions on R+. 


3. Mixtures of Bessel distributions 

To find a distribution vt with the characteristic function of the form (2.8) and 
study the distributional properties of the BAJD, it is inevitable to encounter the 
Bessel distributions and mixtures of Bessel distributions. 

We start with a slight variant of the Bessel distribution defined in [8, p.l5]. 
Suppose that a and /? are positive constants. A probability measure fia,i 3 on 
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l+,S(K+)) is called a Bessel distribution with parameters a and /3 if 

^a,p{dx) = e~°‘6Q{dx) + /— • Ji(2-\/ al3x)dx, (3-1) 


where 6q is the Dirac measure at the origin and Ii is the modified Bessel function 
of the first kind, namely 






{yy 


2^^kl{k + i)r 


r e M. 


Now we consider mixtures of Bessel distributions. Let 7 > 0 be a constant and 
define a probability measure on ]R_|_ as follows: 

/*oo 

ma,p,^{,dx) := / p,at,p{dx)-—e~*dt. 

Jo 1 ( 7 ) 

Similar to [ 8 ] we can easily calculate the characteristic function of pia,p and 

^ck,/ 3 , 7 - 

Lemma 3.1. For u €14 we have: 


POO 

(z) / e'^^Ha,p{dx) = 

Jo 


(ii) / e“^m, 


,P,'ridx) = ( 
Proof, (i) li u €14, then 


a + 1 a + 1 1 — 


-)i 

• u/ 


pOC CO 

=e-“ + e-“ / • V 

00 

=e-“ + a/3e-“ ^ 


dx 


-n^O 


/c =0 




(u-P). , 

A:!(fc + 1)! 
a/3 \fe+i 1 


fc =0 


=e 


P — uJ {k + 1)! 

1 
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(ii) For u vje get 


^0 

= ( 
^ ( 




0 r(7) 

-/3 + (a + l)u\-y 


u — (3/ 
-P + u y 


\-7 ^ / -P + U \ 
) \ —/3 + (a + l)u) 

{{a + l)u- (3) + -^- /3n -7 


)■ 


—j3 + (a + l)u 

/I a 1 \7 

Va + l'^a + 1 1- ^ -u) 


Lemma 3.2. (i) The measure rria^p^-f can be represented as follows: 

ma,l3,^{dx) = So+ ga,p,j{x)dx, x>0, 


where 


y,l3,-fix) 'y ] 


a^r(fc + 7 ) 

^ (a + l)^+'>T( 7 )fc! 


r(a::; k,P), x > 0 , 


(3.2) 


□ 


(3.3) 


(3.4) 


andT{x; k, /3) denotes the density function of the Gamma distribution with param¬ 
eters k and /3. 

(ii) The function ga,/ 3 ,-yix) defined in (3.4) is a continuous function with variables 
{a, (3, 7 , a:) G D := (0, 00 ) x (0, 00 ) x (0, 00 ) x [0, 00 ). 

Proof, (i) We can write 

roc 

rua^pydx) = / tiat,p{dx)--^e~*^dt 

Jo ^ ( 7 ) 

/»00 . / I 17 — 1 

= J ye~°'* 5 o{dx) 3 -/ 3 e~°‘*~^^^ — ■ Il{ 2 ^/cd^)dxj-^^t 


^dt 


(^YSoidx)+ 


dtdx 


1 p 00 


g-(a+l)i^7+fc 

r(7)fc!(fc + i) 


-dt) 


dx 




a'=+ir(A: + 7 + 1) 


fc =0 

00 ,,k 




^ (a + l)^+^r(7)fc! 
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(ii) By the definition of ga,p,-y{x) we have 

Hxt 


9a,P,i{^) — I 

Jo 


Pe~ 


ajJe 


_ ^7-1 

/i( 2 y' atj3x) „, ^ e~*dt 


*dt 


r(7) 

. {at(3x)^ . n 


r(7) 


:e 

fc =0 


k\{k + l)\' 


Suppose that (ao) Po, 70 ) xq) G D and <5 > 0 is small enough such that 70 — 5 > 0, 
ap — ^ > 0 and /3o — > 0. Then for (a, /3, 7 , x) G Ks with 

Ks ■■= {(a,/ 3 , 7 ,x) G D : max{|a - ap], |/3 - /3o|, I 7 - 7 o|, \x - xpl} < (5} 

we get 


:e 


{at/3x)^ 


^g-(a+l)t-/3, 

r(^) v^^fc!(fe + i)!^ ^r(7) 




al3t^ (^+i)t-8x(^ {atf{l3xf 


E 

A;=0 


(fc !)2 


-) 


apt^ rct+i)t-8x at Bx < t 

-r(7) ■ -r(7) 

<C5 + r'>+^e-‘l(i,o,)(t)) (3.5) 

for some constant cs > 0 , since is continuous and thus bounded for 
(a,/3,7,a:) G Ks- If {an, Pn,ln,Xn) -t (ap,/3p, 7 e, xp) as n -)> 00 , then by domi¬ 
nated convergence we get 

lim gan,l3n,ln{^'n) = 5c(o ,/3o , 7 o (^0 ) I 


namely ga,8,'y{x) is a continuous function on D. 


□ 


Remark 3.3. If we write Sq = r(0,/3), namely considering the Dirac measure (5p 
as a degenerated Gamma distribution, then the representation in (3.3) shows that 
the measure TOa,/ 3,7 is a mixture of Gamma distributions r(fc,/3), k G namely 


(t^)>(m)+i: 


Too 


k=l 


a'^T{k + 7 ) 

(a + I)^'+'>T( 7 )fc! 


r(fc,/3). 


4 . Transition density of the BAJD 

In this section we shall derive a closed form expression for the transition density 
of the BAJD. We should mention that in [4, Ghapter 7] the density functions of 
the pricing semigroup associated to the BAJD was derived for some special cases. 
Essentially, the method used in [4] could be used to derive the density functions 
of the BAJD in the case where c/(a — a^d/2) G Z. Here we proceed like [4] but 
deal with more general parameters. In order to do this, we first find, by using the 
results of the previous section, a probability measure vt on 1 R+ whose characteristic 
function satisfies ( 2 . 8 ). 

We recall that the BAJD process X = (Xt)t>o is given by (l.I). We distinguish 
between three cases according to the sign of A := a — cP’dj^. 
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4.1. Case i): A > 0. From (2.5) we know that 

= (1 - ^uil - • exp ( --) 


•( 


d- 


^du 


2a 


T -1)” 


d — u 


)' 


(4.1) 


The product of the first two terms on the right-hand side of (4.1) coincides with 
the characteristic function of the CIR process Z = {Zt)t>o defined in (2.6). It is 
well-known that the transition density function of the CIR process is given by 


f{t,x,y)=pe “ ^ Iq{2{uv)^^) 


(4.2) 


for t > 0, a; > 0 and y > 0, where 

2a 


P =- 


cr^ ^1 — e y 


V =py, 


u = pxe 


2a9 , 

q = - 1, 


and Iq{-) is the modified Bessel function of the first kind of order q, namely 

oo 


fe=0 


k\T{q -I- fc -I- 1) ’ 


We should remark that for a; = 0 the formula of the density function f{t,x,y) 
given in (4.2) is not valid any more. In this case we have 

P ^ — V 


/(CO, 2/) = 


r(<z + i) 


(4.3) 


for t > 0 and ?/ > 0. 
Thus 


f /(Cai,i/)e“^d 2 /= (l - ^u(I - e “*)) ^ • exp 
Now we want to find a probability measure ut with 


(iTi 


xue 


I - f-ir(I - e-“‘) 


)• 




(dy) =( 


d-' 


d 


/ d — uLi{t) \ 

V d — u ) 

-I- (l - Li{t)) ^ 


(4.4) 


where Li{t) := exp(—at) -I- (T^d(l — exp(—at))/(2a). If such a measure Ut exists, 
then the law of Xt can be written as the convolution of the law of Zt and Uj. 

Comparing the characteristic functions (3.2) and (4.4), it is easy to see that we 
can seek the measure vt as a mixture of Bessel distributions. More precisely, we 
define 




(4.5) 
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with 

< := (4.6) 


Then the characteristic function of uj coincides with (4.4). Since the probability 
measure is of the form (3.3), it follows now from (2.7), (4.1) and 

(4.4) that the law of Xt is absolutely continuous with respect to the Lesbegue 
measure and its density function p(t, x, y) is given by 

/ 1 \7l 

p{t,x,y) = j f{t,x,y) + J f{t,x,y - (4-7) 

for t > 0, X > 0 and y > 0, where the function g is defined in (3.4). 


4.2. Case ii): A < 0. Similar to the case (i), it suffices to find a probability 
measure ut with 


f e"Vt(dy) = (- 

Jr+ ^ 

- ( 


d — u 
d — u 


.)7 


ct '^ d 


2a 

/ d — U N “ 

\d — (t)u) 


ue 






d — Li{t)u 
1 


Li{t) 


(1-^) 




Since A = a — cr^d /2 < 0, therefore a^dj^a > 1 and 

Ti(t) = e-“‘ + ^-(l-e-“*) >1. 
According to the formula (3.1), we can choose 

Ut ^C(2(t),/32(t),72 

with the parameters a 2 , P 2 and 72 dehned by 

a2{t) := Li{t) - 1 
/32 :=d 
72 := 




(4.8) 


(4.9) 


Similar to the case (t), the transition densities p{t,x,y) of X is given by 
/ 1 \ 72 py 

p{t,x,y) = [— - —) f{t,x,y)+ f{t,x,y - z)ga^^t),p 2 ,i 2 iz)dz (4.10) 

\l + a2[t)y Jq 

for t > 0, a: > 0 and y > 0, where the function g is defined in (3.4). 
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4.3. Case iii): A = 0. In this case we need to find a probability measure vt 
with 

'cu{l — e““‘)' 


e'^yvt{dy) = exp 


a{d — u) 


According to the formula (3.3) we can take vt as a Bessel distribution yLa 3 {t),p^ 
with the parameters 03 (t) and (3^ defined by 


iasit) :=f(l-e-“‘) 

\/33 := d. 

Thus in this case the transition densities p{t, x, y) of X is given by 

p(t,x,y) = j f{t,x,y - ^^^^^Ii{ 2 y/a 3 {t)P 3 z)dz 

+e~°‘^^*^fit,x,y) 
for t > 0 , X > 0 and y >0- 

Summarizing the above three cases we get the following theorem. 


(4.11) 


(4.12) 


Theorem 4.1. Let X = {Xt)t>o be the BAJD defined in (1.1). Then the law of 
Xt given that Xq = x > 0 is absolutely continuous with respect to the Lesbegue 
measure and thus posseses a density function p{t,x,y), namely 

Px{Xt e A) = [ p{t,x,y)dy, t>0, A e 6(]R+) 

J A 

According to the sign of A = a — a‘^d/2, the density p{t,x,y) is given by (4-7), 
(4-10) and (4-12) respectively. 

Although the density functions in (4.7), (4.10) and (4.12) are essentially differ¬ 
ent, they do share some similarities. In the following corollary we give a unified 
representation of p{t, x, y). 

Corollary 4.2. Irrelevant of the the sign of A = a—a^d/2, the transition densities 
p{t, X, y) of X can be expressed in a unified form as 

p{t,x,y) = L{t)f{t,x,y) + [ f{t,x,y-z)h{t,z)dz, (4.13) 

Jo 

where L{t) is continuous function in t > 0 which satisfies 0 < L(t) < 1 for t > 0, 
the function hft, z) is non-negative and continuous in (t, z) G (0, 00 ) x [0, 00 ) and 
satisfies h(t, z)dz = 1 — Lft). 


5. Positive Harris recurrence of the BAJD 

It was shown in [2] (see also [13]) that the semigroup of any stochastically con¬ 
tinuous affine process on the canonical state space R™ x R” is a Feller semigroup. 
Define the semigroup of the BAJD by 

Ttfix):=[ p{t,x,y)f{y)dy, 


(5.1) 
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where / : R+ —R is bounded. Since the BAJD process X is stochastically 
continuous and afhne, thus {Tt)t>o is a Feller semigroup. 

To show the positive Harris recurrence, we need hrst to prove the regularity 
property of BAJD. To this aim, we first analyse the continuity properties of the 
integral which appears on the right hand side of (4.13). 

Lemma 5.1. Let f{t,x,y) be the transition density of the CIR proeess given in 
(4-2) and h{t,z) be the same as in (4.13). Then the funetion F{t,x,y) defined by 

F{t,x,y):= [ f{t,x,y-z)h{t,z)dz (5.2) 

Jo 

is continuous with variables if, x, y) G (0, oo) x [0, oo) x [0, oo). Moreover if M > 1 
is a constant, then 

\Fit,x,y)\<Cy"i^ (5.3) 

for all 

{t, x, y) G Km ■■= {{t, x,y) : ^ <t < M, 0<x<M, 0<y<^}, 
where C > Q is a constant which depends on M. 

Proof. For simplicity we set q := 2a9fa'^ — 1 as in (4.2). Since h{t, z) is continuous 
in {t,z) G (0,oo) X [0, oo), thus there exists a constant ci > 0 depending on M 
such that 

\h{t, z)\ < Cl for ^<t<M, 0<z<^. (5.4) 

Therefore if {t,x,y) G Km, we have 

\F{t,x,y)\<Ci [ f{t,x,y - z)dz. (5.5) 

Jo 

According to (4.2) and (4.3) we have 

\f{t,x,y - z)\ < C 2 \y - z\'^ if {t,x,y - z) € Km and y ^ z, (5.6) 

where C 2 > 0 is a constant depending on M. 

It follows from (5.5) and (5.6) that 

\F{t,x,y)\<ciC 2 [ \y - z\‘^dz =< C 3 y‘>"^^ (5.7) 

Jo 9 + 1 

for {t,x,y) G Km, if we set C 3 := ciC 2 f{q + 1). Thus (5.3) is proved. Noting that 
F{t,x,0) = 0 for t > 0 and a; > 0, the continuity of the function F at points 
(t,x,0) is an immediate consequence of the estimate (5.7). 

We now proceed to prove the continuity of F at other points. Suppose that fp > Oj 
aiQ > 0 and yo > 0 are fixed. Let e > 0 be arbitrary. We choose > 0 small 
enough such that j/q — 25i > 0 and fp — Ji > 0. As in (5.4) and (5.6) there exist 
constants 04,05 > 0 , which depend on Ji, such that 

\h{t,z)\<Ci for t € [to - Si,to + Si], ze[0,yo + Si] (5.8) 


and 


\f{t,x,y- z)\ < C5\y - z\‘^ 


( 5 . 9 ) 
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for i G [io — <5i, to + <^i], x G [0, xo + (5i] and 0<i/ — 2:<?/o+<5i- Set 

Ks 2 ■= [to — <52,to + < 52 ] X [0,a:o + <52] x [?/o — <52, Vo + <52]- 

We choose <52 > 0 small enough such that S 2 < <5i and 0405 ( 362 )’^^^/{q + 1) < e/3. 
If (t, X, y) G then it holds 


f{t,x,y- z)h{t,z)dz 


yo-262 


< C4C5 


iVa-252 


{y- zYdz 


= ^^(y-yo + 2,52)®+^ 

< ^(3<52)«+' < (5.10) 

g + 1 3 

If (t, X, y) G and ^3 < z <yo — 282 , then 82 < y — z < yo + 82 and by (5.8) and 
(5.9) we have 


\f{t,x,y - z)h{t,z)\ < C 4 C 5 |y- 2 :|'’ 

< C4C5(|<52r + |?/o + <52r). (5.11) 


Since for fixed z G [0, j/o ~ 2 ^ 2 ] the function f{t,x,y — z)h(t,z) is continuous in 
(t,x,y) G Ks 2 , it follows from (5.11) and dominated convergence theorem that 

pyo—252 

F 2 {t,x,y):= f{t,x,y-z)h{t,z)dz 

Jo 

is a continuous function in (t, x, y) G Ks 2 - This implies the existence of a constant 
8 with 0 < 8 < 82 such that 

\F2{t,x,y) - F2{to,xo,yo)\ < |, (5.12) 

if (t, x, y) G Kg := [to - 5, to + (5] x [0 V (xo - <5), a:o + 5] x [yo - i5, yo + <5]. Thus it 
follows from (5.10) and (5.12) that 

\F{t,x,y) - F(to,a:o,yo)l < e 


for (t,x,y) G Kg. The continuity of the function F at (to,Xo,yo) is proved. □ 


Following [10, Chapter 20] we give the definition of a regular Feller process on 
R+. 

Definition 5.2. Consider a Feller process Y = {Yt)t>o with state space (R+, S(IR+))| 
and distributions Px, x G K+. The process is said to be regular if there exist a 
locally finite measure p on 1 R+ and a continuous function {t,x,y) i-A pt{x,y) > 0 
on (0, 00 ) X such that 

Px{Yt G B} = I pt{x,y)p{dy), x G M+, BgB{R+), t > 0. 

Jb 

The measure p is called the ’’supporting measure” of the process. It is unique up 
to an equivalence (see [10, page 399]). 


Proceeding as in the proof of [9, 
on (]R 4 _, S(]R+)) as 


Proposition 2.2], we define a new measure 77 


( 5 . 13 ) 


rj{dx) := K{x)dx, 






POSITIVE HARRIS RECURRENCE AND EXPONENTIAL ERGODICITY OF THE BAJD 13 


where 



0 < a; < 1, 

X > 1. 


Then the transition densities of the BAJD process with respect to the new measure 
r] is given by 


P{t,x,y) 


P{t,x,y) 

K(y) 


t > 0, X >0, y > 0. 


(5.14) 


Theorem 5.3. The transition densities pit, x, y) of the BAJD process with respect 
to the measure rj satisfies 


0 < p{t, X, y) < oo, t > 0, X > 0, y > 0 

and is continuous in {t,x,y) € (0, oo) x [0, oo) x [0,oo). Thus the BAJD is a 
regular Feller process on the state space R+. 


Proof. From (4.13) we can write the transition density of the BAJD process with 
respect to the measure rj as 


p(,t,x,y) 


m 


f{t,x,y) 

«(y) 


L{t)f{t,x,y) 


F{t,x,y) 

K(y) 

F{t,x,y) 


(5.15) 


where f{t,x,y) is the transition density of CIR with respect the new measure ry 
and F is defined in (5.2). We have already shown in [9, Proposition 2.2] that 

0 < f{t,x,y) < oo, t > 0, X > 0, y > 0, (5.16) 


and {t,x,y) i—>■ f(t,x,y) is a continuous function on (0,oo) x 

From Lemma 5.1 we know that the function F(t,x,y) appearing in the second 

summand in (5.15) is continuous on (0,oo) x and 


\F{t,x,y)\ < Cy^, if {t,x,y) € Km, 


where C > 0 is a constant depending on M. Now it is clear that 

0 < hm F ^ ^11 ^ 

y^O K{y) y^O 

Since F{t,x,0) = 0 and L{t) is continuous in t > 0, it follows that the function 
p{t,x,y) is continuous at points (to,a;o,0), if to > 0 and Xg > 0. The continuity 
of the function p(t, x, y) at other points is also clear, because all the functions 
appearing in (5.15) are continuous and 0 < K{y) < oo for y > 0. Noting that 
0 < L{t) < 1 for t > 0, we get 

0 < p{t, x,y) < oo for all (t, x, y) S (0, oo) x R^. 

Therefore the BAJD process is a regular Feller process with 77 as the supporting 
measure. □ 


For the following definitions we refer to [17] . 
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Definition 5.4. Consider a time-homogeneous Markove process Y = (Yt)t>o with 
the state space K+ and distributions Px, x G K+. 

(i) Y is said to be Harris recurrent if for some cr-finite measure /i 


Px{J lA{Ys)ds = = 1, (5.17) 

for any x G K+ and A G H(R+) with /^(A) > 0. It was shown in [7] that if Y 
is Harris recurrent then it possesses a unique (up to a renormalization) invariant 
measure. If the invariant measure is finite, then the process Y is called positive 
Harris recurrent. 

(ii) Y is said to be uniformly transient if 


supifx 

X 


lK(Ys)ds 


< oo 


(5.18) 


for every compact K C R+. 


Lemma 5.5. The BAJD is not uniformly transient. 

Proof. Let m > 0, K := [0, m] and x G (0, oo) be fixed. Then 


Px 


‘-xo 


^[O.m] (.^t)dt 


/ Px\)-[0,m]{Xt)\dt 

h 

poo pm 

/ / p{t, X, y)dydt 

/o ^0 

pm poo 

/ dy p{t,x,y)dt. 
fo Jo 


From (4.13) and (5.2) we know that 

p{t, X, y) = L{t)f{t, x, y) + Fit, x, y). 
It follows from (4.7), (4.10) and (4.12) that 


0 < lim L(t) = I (‘2a‘^)“ ^ ’ if A 7 ^ 0 , 

[e-^, ifA = 0. 

Since Fit,x,y) is non-negative, thus there exists large enough T > 0 such that 
p{t, X, y) > Xf{t, X, y) for t > T, 

where A > 0 is a constant. Let e > 0 be small enough. According to [9, Lemma 
2.4] we know that for any a: > 0 and y G [e, m] it holds 


pOO 

/ f{t,x,y)dt = oo. 
Jo 


Therefore 


pm poo 


Fx / l[o,m](At)dt = / dy p{t,x,y)dt 
^Jo ^ Jo Jo 

pm poo 

> J dy J Xf{t,x,y)dt 

This proves that the BAJD is not uniformly transient. 


□ 
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Theorem 5.6. The BAJD is Harris recurrent. 

Proof. We have shown that the BAJD is a regular Feller process with the measure 
77 , which is defined in (5.13), as a supporting measure. It follows from Lemma 5.5 
and [10, Theorem 20.17] that the BAJD is Harris recurrent and the measure 77 
satisfies (5.17). □ 


Remark 5.7. Since the BAJD is a Harris recurrent Feller process, it follows from 
[10, Theorem 20.18] that the BAJD has a locally finite invariant measure tt which 
is equivalent to the supporting measure 77 . Moreover every a -finite invariant 
measure of the BAJD agrees with tt up to a renormalization. The existence and 
uniqueness of an invariant probability measure for the BAJD has already been 
proved in [11] (see also [14]). Thus we can assume tt to be a probability measure. 
The characteristic function of tt was given in [14] and has the form 



(1 

(1 



if A ^ 0, 


if A = 0. 


(5.20) 


Corollary 5.8. The BAJD is positive Harris recurrent. Its unique invariant 
probability measure tt is absolute continuous with respect to the Lehesgue measure 
and thus has a density function l{-), namely T^{dy) = l{y)dy, y € M-f.. // A 7 ^ 0, 
then we have 




2a9 

cr^ ’ 2 a 



2a0 

’ ’ 2 a 


)h{z)dz, 


2/ > 0, 


where T{jj]2a9j{2a)^ denotes the density function of the Gamma distribu¬ 
tion with parameters 2a9/a'^ and jifla'), and 


h{z) 




if A > 0, 
if A < 0, 


with g defined in (3.4). 4/A = 0, then we have 


l{y) = e -r(y; 


2a9 

(T^ ’ 2 a 


ry 

)+ r(7/ 

Jo 


- Z] ■ 


2a9 (7^ 


)d( 


— — —dz 


C 

adz 



dz 


for y >0. 


Proof. We apply the same method which we used in Section 4 to find the transition 
densities of the BAJD. Since the characteristic function of tt is given by (5.20) and 
noting that the first term on the right hand side of (5.20) corresponds to the 
characteristic function of a Gamma distribution, we can represent the measure 
TT as a convolution of a Gamma distribution with a probability measure v. If 
A = 0 the measure u is a Bessel distribution, otherwise it is a mixture of Bessel 
distributions. By identifying the parameters of the Gamma distribution and the 
Bessel or mixture of Bessel distributions, we get an explicit formula of the density 
function of tt. □ 
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Corollary 5.9. Let X = {Xt)t>o be the BAJD defined by (1.1). Then for any 
f G Sb(K.+) we have 


I f f(Xs)ds f f{x)TT{dx) a.s. 
t Jo JRi 


as t ^ OQ, where tt is the unique invariant probability measure of the BAJD. 

Proof. The above convergence follows from Corollary 5.8 and [10, Theorem 20.21]. 

□ 


6 . Exponential ergodicity of the BAJD 

Let II • \\tv denote the total-variation norm for signed measures on namely 

M\tv= sup {|/j(A)|}. 

AgB(R+) 

The total variation norm for signed measures on ]R_|_ is a special case of the the 
norm || ■ ||b, which is defined by 



for a function h on R+ with h > 1. Obviously it holds ||/j||Ty < IImIUi given that 
h>l. 

Let P*{x, •) := Px{Xt £ •) be the distribution of the BAJD process X at time 
t given that Xq = x with a; > 0. In this section we will find a function h>l such 
that the BAJD is /i-exponentially ergodic, namely there exist a constant fi G (0,1) 
and a finite-valued function B{-) such that 

\\P\x,-) - Tr\\h < B{x)l3\ Vt > 0, a; G M+, 

where tt is the unique invariant probability measure of the BAJD. 

We first show the existence of a Foster-Lyapunov function, which is a sufficient 
condition for the exponential ergodicity to hold. Let A denote the generator of 
the BAJD, then 

■^9{x) = ^<T'^xg"{x) + {ad - ax)g'{x) + cd [ {g{x + y) - g{x))e~‘^ydy 
for g G D(A). 

Lemma 6.1. The function V{x) = exp( 7 a:) with small enough j > 0 is a Foster- 
Lyapunov function for the BAJD, namely there exist constants k,M € (0,oo) such 
that 


AV(x) < —kV (x) M, \/x € IR+ 




POSITIVE HARRIS RECURRENCE AND EXPONENTIAL ERGODICITY OF THE BAJD 17 

Proof. We have 

AV{x) + {aO - ax)-fe^^ + cd [ 

=-cr^7^xe'^'^ + (aO — ax)^e^^ + 

2 d — 7 

=76'’''“ • + ad + ^ . 

If 7 > 0 is small enough, then (7^7/2 — a < 0 and there exists xo > 0 with 

/ 1 2 \ /I ^ ^ 

[-a 7 - a)xo + + - -=-. 

^2 ^ a — 7 7 

Thus we have for x G [xq, oo) 

AV{x) = 76'’'“ • f(^cr^7 — a)x + a 9 + \ < —e'’'® 

and for x £ [0,a:o] 

AV{x) = 76'’"“ • f (-a^7 — a)a; + ad + - 7 —) — 76”’’'^“ ' (ad + - 7 —)■ 

V '2 ^ a — 7/ a — 7 

It follows for all x £ R+ 

AV{x) < -e<^ + 6 ^"““ + 76^'“° • (ad + - 7 —) < -^(2’) + ^ 

a — 7 

with 

M := 6 ^'““ + 76^'““ • (ad + -^) < 00. 


Lemma 6.2. Let t/ie constants 7 , /c and M, as well as V{x) 
same as in Lemma 6.1. Then the BAJD satisfies 


E,[V{Xt)] < e-^W{x) + ^ 


or equivalently 


f V (y)p(t, X, y)dy < e '"*V{x) + 
Jr , 


M 

~k 


for all X £ M+, t > 0. 


exp ( 7 a;) be the 


Proof. Let V(x) = exp( 7 x) and g(x,t) '■= r(x) • exp(fct) = exp( 7 x + kt), where 
the constants 7 and k are the same as in Lemma 6.1. Then gx = 'yexjp{'jx + kt), 
gxx = 7 ^ exp( 7 a; + kt) and gt = kexp{jx + kt). By applying Ito’s formula and 
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then taking the expectation, we get for all x gk+, t >0, 


e'^*E,[V{Xt)]-V{x) 

=E,[g{Xt,t)]-EMXoM 


<E, 

=Ex 



■AV{Xs) + ke^‘’-V{Xs))ds 


^0 


Xle'^^ds 


kV{Xs) +M)+ ke^^ ■ V{Xs))ds 


Thus for X G M+, < > 0, 


E:,[V{Xt)\ < e-’^*V{x) + 


□ 


Applying the main results of [18] and Lemma 6.2, we get the following theorem. 

Theorem 6.3. Let h{x) := 1 + exp( 7 a;) with the eonstant 7 > 0 small enough. 
Then the BAJD is h-exponentially ergodic, namely there exist eonstants P G (0,1) 
and C G (0,oo) sueh that 

<C(eT'" + l)/3*, t>0, a: G R+. (6.1) 

Proof. Basically, we follow the proof of [18, Theorem 6.1]. For any (5 > 0 we 
consider the i5-skeleton chain := Xns, n G Z_|_. Then (F(f)nez+ is a Markov 
chain with transition kernel p(5, x, y) on the state space 1 R+ and the law of the 
Yn (started from Iq = a;) is given by P^"(a;,-). It is easy to see that invariant 
measures for the BAJD process {Xt)t>o are also invariant measures for (F)f)ngz_,.. 
Thus the probability measure tt in Corollary 5.8 is also an invariant probability 
measure for the chain (F)f)raez+- 

Let V{x) = exp( 7 a;) be the same as in Lemma 6.2. It follows from the Markov 
property and Lemma 6.2 that 

A,[M(y„+i)lro,yi,--- ,y„] = [ Viy)p{S,Y^,y)dy < e-^^V{Y^) + 

where k and M are positive constants. If we set Vq := V and I 4 , := V(Y„), n G N, 
then 

E4Vi] < e-^^Voix) + ^ 

and 

E4Vn+i\Yo,Yir-- ,Y^]<e-^'^Vn + ^, n G N. 

Now we proceed to show that the chain (L)f)nez+ is A-irreducible, strong ape¬ 
riodic, and all compact subsets of K+ are petite for the chain (F(f)nez+- 
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“A-irreducibility”: We show that the Lebesgue measure A on R+ is an irre- 
ducibility measure for (F^)„gz+- Let A G ;B(R+) and A(^) > 0, then 


P[Yi G A\Yo = x]= Px[Xs G A] = [ p{S,x,y)dy > 0, 

Ja 

since p{S, x,y) > 0 for any x G R+ and y > 0. This shows that the chain (F^)nGZ+ 
is irreducible with A being an irreducibility measure. 

“Strong aperiodicity” (see [16, p.561] for a definition): To show the strong 
aperiodicity of we need to find a set C G S(R+), a probability measure 

u with u(C') = 1, and e > 0 such that 


and 

L{x,C) > 0, 

X G M+ 

(6.2) 


PxiYi G A) > e • u(A), 

X G C, A G 6(M+), 

(6.3) 

where L{x, C) 

:= Vx{Ym G C for some m 

G N). To this end set C 

:= [0,1] and 


g{y) := inf^-gjo i] p(i5, a:, y), y > 0. Since for fixed y > 0 the function p{5,x,y) 
strictly positive and continuous in x G [0,1], thus we have g{y) > 0 and 0 < 
I(o,i]9{y)dy < 1. Define 


'(A) := 


9{y)dy, A G B{R+). 


1(0,1] 9iy)dy JAnio,!]' 

Then for any x G [0,1] and A G S(R+) we get 

Pa:iYiGA)= / p{5,x,y)dy> / g{y)dy = u{A) / y(y)dy, 
JA JAnfo.il aio.il 


so (6.3) holds with e := g{y)dy. 
Obviously 


L(x,[0,l])>P,(yiG[0,l]) = P,(X 5 G[ 0 ,l])= [ p{d,x,y)dy >0 
for all X G 1R+, which verifies (6.2). 

“Compact subsets are petite”: We have shown that A is an irreducibility mea¬ 
sure for (y^)„gz_^. According to [16, Theorem 3.4(ii)], to show that all compact 
sets are petit, it suffices to prove the Feller property of (F)f)„ez+, but this follows 
from the fact that (1 ^)„gz+ is a skeleton chain of the BAJD process {Xt)t>Q, 
which possess the Feller property. 

According to [16, Theorem 6.3] (see also the proof of [16, Theorem 6.1]), 
the probability measure tt is the only invariant probability measure of the chain 
(l)f)„gz+ and there exist constants /3 G (0,1) and C G (0,oo) such that 

||P^"(x,-)-T|U<C'(eT'" + l)/3", nGZ+, x G R+, 
where h := V + I = exp( 7 x) + 1. 

Then we can proceed as in [18, p.536] and get the inequality (6.1). □ 


Since ||/i||Ty < IImIUj it follows immediatly the follwing corollary. 
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Corollary 6.4. The BAJD is exponentially ergodic, namely there exist constants 
/? G (0,oo) and C G (0,oo) such that 

\\P\x,-)-tt\\tv <C{e^^+ 1)I3\ Vt > 0, x G K+. (6.4) 

Remark 6.5. In [9] we proved that the CIR process is positive Harris recurrent. If 
we allow the parameter c = 0, then all results of this section still hold and thus are 
also true for the CIR process. In particular Theorem 6.3 is also true for the CIR 
process. In this case the unique invariant probability measure of the CIR process 
is the Gamma distribution r(2a0/(T^, cr^/(2a)) and has the characteristic function 
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